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The decay width of the process n — » 7r°77 is calculated in the framework of the Nambu-Jona- 
Lasinio model. The momentum dependence of quark loops is taken into account. Three types 
of diagrams are considered: quark box, scalar(ao) and vector(p, w) pole diagrams. The obtained 
estimations are in satisfactory agreement with recent experimental data. 



o 



The investigations of the process rj — ► 7r°77 have a long history. The experimental studies of this process began in 
1966 1 |2|. The first experimental results led to a large value of the branching ratio of the process. The theoretical 
estimates obtained in the vector dominance model (VDM) jij, nonlinear chiral theory 2 Q and lately in the chiral 

■ quark model IE predicted noticeably lower value. 

A real breakthrough in the investigation of this process happened in the experiment GAMS in 1981 at Protvino 

■ Q where the large energies of the produced 77-mesons dramatically suppressed the background. During a subsequent 
reanalysis, the value r, ; ^^ 77 = 0.84 ± 0.18 eV was obtained 3 j9j . Lately, the SND collaboration in the experiment 
VEPP-2M confirmed this value to be an upper limit 1 eV for the width of the process 0. In 2005, the results obtained 
by the Crystal Ball collaboration at BNL AGS were published; these results r^_» 7r77 = 0.45±0.12 El were noticeably 
smaller than those reported by the GAMS collaboration. 

From a theoretical point of view this process was investigated in many theoretical models: VDM model Q , nonlinear 
chiral theory jjj L di fferent quark models 0, El El El j resonance exchange models EIE|, the chiral perturbation 
__, i theory (ChPT) | la . I Hi ll 71 Il8l Il9|| . and chiral unitary approach 20] . In ChPT, the main contribution comes from the 
I terms of the order 0(p 6 ) of low energy expansion because the tree terms of the order 0(p 2 ) and 0(p 4 ) are absent and 
1 one-loop contributions of the order 0(p 4 ) are very small. The counterterms of the order 0(p 6 ) are not determined from 
the theory itself and should be fixed using experimental information, from the assumption of meson saturation (vector 
meson exchange giving the dominant contribution) or calculated from the model (NJL for example) . In [Tol , the meson 
saturation approach was adopted, which gave T l; ^ 7r77 = 0.18 eV; too small, compared to the experimental value. But, 
\^ \ keeping the momentum dependence in the vector meson propagators gives an "all-order" estimate of about 0.31 eV 
|l5j |. in agreement with the old VDM prediction 0. Taking into account the scalar and tensor meson contributions 
(the signs of which cannot be unambiguously determined within this approach) and the one-loop contribution at 
0(p 8 ), the final estimate of El i s T^^^v = 0.42 ± 0.20 eV, in a satisfactory agreement with the recent Crystal 
Ball result. This result is confirmed in 16], where the 0(p 6 ) counterterms are calculated in the framework of the 
Qh. NJL model with the result 0.58 ± 0.3 eV. However, the same counterterms obtained from the NJL model by different 
methods lead to 0.1 eV [H| and 0.27±g;o? eV [p|. 

The "all-order" estimations in El are a signal that the preservation of full momentum dependence is highly 
desirable. Note that in 0, the simple NJL model is used without taking into account the momentum dependence of 
quark loops. Then, in a quark models pdlll2|. the full momentum dependence of the quark box diagram is considered 
whereas the diagram with the intermediate scalar ao(980) is ignored. The vector sector of the model has not been 
taken into account as well. 

In the present work, the process rj — > ^77 is calculated in the framework of the NJL model with scalar-pseudoscalar 
and vector-axial-vector sectors. The contribution of the quark box loop is considered together with the contributions 
of the diagrams with scalar and vector intermediate mesons (as in 0,0)- The momentum dependence of the quark 
loops and pseudoscalar-axial- vector transitions are taken into account, following pll I22I |23| . 
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1 Excellent review of theoretical and experimental works can be found in 

2 Note that similar result for the width of the order 10 — 2 eV lately obtained in ChPT for pion-loop contribution at the level 0(p 4 ). 

3 Notice that this result is consistent with those obtained in the NJL model 0\. 
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I. THE (7(3) x 1/(3) NJL MODEL 



The t/(3) x U(3) NJL model with scalar-pseudoscalar and vector-axial-vector sectors is used in the present work. 
To solve the J7a(1) problem, the six-quark t'Hooft interaction is added to the Lagrangian of the model p3.l25| 



G 



C = q(id - m°)q + — ^2[(q\ t q) 2 + (qij5 hqf] + ~7T ^[(Qln^iOf + {qi^l^Kq) 



where \ (i=l, 



-K (det[f (1 + 7s)g] + det[g(l - 75)9]) , 



are the Gell-Mann matrices and A 



(1) 



1, with 1 being the unit matrix; m° is the current 



quark mass matrix with diagonal elements mj, rnP s (m° ~ m d)i G and Gy are the scalar-pseudoscalar and 
vector-axial-vector four-quark coupling constants; K is the six-quark coupling constants. The six-quark interaction 
can be reduced to an effective four-fermion vertex after the contraction of one of the quark pairs. The details are 
given in appendix A. 

Light current quarks transform to massive constituent quarks as a result of spontaneous chiral symmetry breaking. 
Constituent quark masses can be found from the Dyson-Schwinger equation for the quark propagators (gap equations) 



to,, = TO,, 



8m u Gh(m u ) + 32m u m s K Ii(m u )Ii(m s ) 



m s = + 8m s GIi(m s ) + 32K (m u Ii(m u )Y 



(2) 



where I\{m) is the quadratically divergent integral. The modified Pauli-Villars (PV) regularization with two sub- 
stractions with same A is used for the regularization of divergent integrals 4 (see [2l], I22I l23l H(|). In this case the 
quadratically and logarithmically divergent integrals I\ (to) and 1% (to) have the same form as in the four- momentum 
cut-off scheme 



him) = J-; 



A 2 

A^ - m z In ( — + 1 

TO 



T( \ Nc 



m4 

Ul 



1 - 



2x - 1 

TO 

1 + A 2 " 



Moreover, the Pauli-Villars regularization is suitable for the description of the vector sector because it preserves 
gauge invariance. 

Masses and vertex functions of the mesons can be found from the Bethe-Salpeter equation. The expression for the 
quark-antiquark scattering matrix is 



f = G + GU(p 2 )f 



1 



g 1 - n( P 2 ) 



(3) 



where G and Tl(p 2 ) are the corresponding matrices of the four-quark coupling constant and polarization loops. The 
particle mass can be found from the equation det(G" 1 — II(M 2 )) = and near the poles the corresponding part of 
the T matrix can be expressed in the form 



T = 



V ®V 



P 



M 2 



(4) 



where V and M are the vertex function and mass of the meson, and V = 7 1^7°. Details of calculations for different 
channels are presented in appendices B, C. Here we discuss only general properties. 

The most simple situation takes place for the vector and the isovector scalar meson with equal quark masses 
(say p and ao). In this case, the coupling constant and polarization operator are just numbers (not matrices). For 
pseudoscalar mesons, additional axial- vector components appear in the vertex function due to the pseudoscalar-axial- 
vector mixing (in the scalar case this transition loop is proportional to the difference of quark masses). An additional 
complication takes place for 77 and rj due to the singlet-octet mixing (or mixing of strange and non-strange quarks 
due to the t'Hooft interaction). Therefore, the vertex function of this meson has four components: strange and 
non-strange pseudoscalar and axial- vector. 



Any function /(m 2 ) of mass m 2 is regularized by using the rule 



/(m 2 ) -» f(m 2 ) - f(rn 2 + A 2 ) + A 2 /'(m 2 + A 2 ). 
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FIG. 1: Diagrams contributing to the amplitude of the process r\ — » 7r 77. 



II. FIXING MODEL PARAMETERS 

The model has six parameters: the coupling constants G, Gv, K, PV cut-off A, and constituent quark masses m u 
and m s . We use two paramctrization schemes. In the first one, the model parameters are defined using masses of 
the pion, kaon, p and r\ mesons and the weak pion decay constant f„. Note that the number of input parameters is 
greater than the number of physical observables by one. This allows us, following 22], to take the mass of the u quark 
slightly larger than the half of the p-meson mass. As a result, we have the following set (set I) of model parameters 

m u = 390 MeV, m s = 496GeV, G = 6.62 GeV~ 2 , G v = -11.29 GeV~ 2 , K = 123 GeV~ 5 , A = 1 GeV. (5) 

The values of the current quark masses m^m® are defined from the gap equations © m° = 3.9 MeV and m° s = 70 
K/m° = 18). 

For this set of model parameters, the two-photon decay width of the 77 meson r,,^ 77 = 0.37 KeV, is smaller than 
the experimental one: T^ 77 = 0.510 ± 0.026 0. 

In the set II the model parameters are fixed in order to reproduce the two-photon decay width of the 77 meson 
instead of its mass (the 77 meson mass in this case M n — 530 MeV) 

m u = 390 MeV, m s = 506 GeV, G = 8.04 GeV" 2 , G v = -11.29 GcV~ 2 , K = 77 GeV" 5 , A = 1 GeV. (6) 
The current quark masses are m° = 3.9 MeV and m° s = 78 MeV (m°/m° = 20). 



III. DECAY 77 -> tt°77 

The general form of the 77 — > n°jj decay amplitude contains two independent tensor structures [26 
T = T^eiei, T"" = A(x u x 2 )(q^ ~ Qi ■ 81 «T) + B( Xl ,x 2 ) 



-(7) 



where p, qi, q 2 are the momentum of the 77 meson and photons, and e 2 are the polarization vectors of the photons, 
and Xi = p - qi/M?. 



The 77 — ► 7r°77 decay width has the form 

(1-1/) /2 xl 



M 5 
256tt 2 



dxo 



A{xi,x 2 ) + -B(x l7 x 2 ) 



[2(xi + ar 2 )+y-l] 



+ i |fl(xi,a;2)| 2 i^ix 2 - [2{ Xl +x 2 )+y- l]] 2 

xf n - (1 - 2xx - y)/2, 4 nax = (1 - 2n - j/)/2(l - 2xx), y = M*/M*. 



(8) 



In the NJL model the amplitude for the 77 — > 7r°77 decay process is described by three types of diagrams (see Fig. 
^): the quark box and exchange of scalar(ao) and vector (p, to) resonances. Let us consider theses contributions in 
detail. 

The scalar meson exchange has the simplest form. It gives a contribution only to A(x\,x 2 ). This contribution 
consists of three parts and can be written in the form (see appendices B and C for the definition of polarization loops 
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TABLE I: rj — > 7r°77 decay width. 



Contribution 


model 1 


model 2 


vector mesons 


0.17 


0.20 


scalar meson 


0.03 


0.12 


vector+scalar mesons 


0.10 


0.12 


box 


0.28 


0.35 


box+vector 


0.78 


0.95 


total 


0.53 


0.45 



and vertex functions): 

A(X!,X 2 ) 



flu 



r(p 2 



ya r]TT (2<7l • <7 2 )5ao77( 2 <?l ' Q2) 
Ga ( l - n™ (2 9l ■ q 2 ) 
1 1-xi 

dxi J dx2 
11 



M 2 



1 

2^ 



qi ■ 92 

m„(l - Axix 2 ) 



X\ + x 2 



Ml 

2 



(p 2 XlCC 2 



A 2 ) 2 (p 2 xix 2 — m 2 ) 



(9) 



9aQ7]7T (P 



-i2N r N 



d\k 



(27T. 



jTr_D {Vao-SUfc + (?i)T/ 7r S' u (fc)K,5 , u (fc - 172)} ■ 



here Tr^ is the trace over Dirac indices, index A in the measure of integration means PV regularization of the integral 
and Sj{p) = (p — rrij)^ 1 . 

The amplitude with the vector meson (p, u>) exchanges consists of two quark triangles of anomalous type (see 
appendix D) and the vector meson propagator. It gives the following contributions 



B(xi,x 2 ) = 
A(xi,x 2 ) = 



E E 

j—p.uj i— 1,2 



ff?m (M 2 ,M 2 (l - 2a; 4 ),0) 5 „ 7 (M 2 ,M 2 (l - 2^,0) 



G 



J2 g ^ (M ^ M ^ (1 

j—p.to i— 1,2 



^-II™ (M 2 (l-2x,)) 

2^, 0) g ^ 7 (M 2 , M 2 (l - 2xj),0)Atf(l - Xj 



(10) 



n™ (Af 2 (i-2^)) 



The box diagram is of a more complicated structure. It consists of three types of boxes (plus three crossed) and 
contains the diagrams with pseudoscalar and axial-vector components of the 7r and rj mesons 



d\k 
(2tt)4 



D 



V 7r S(k)V v S{k+p~ qi - q 2 )^ u S(k +p- qi)^ ll S{k + p) + 
+V„S(k)V n S(k + q 2 ) lu S{k +p- qi )^S(k + p) 

+V n S(k)j„S(k + q 2 )-f ll S{k + qi+ q 2 )V v S(k+p) + {q x <-> q 2 ,fi <-> v} 



(11) 



We calculate these diagrams numerically. In order to check the integration procedure, we calculate all coefficients of 
different tensor structures and verify if they have gauge invariant form JJJ. 

The obtained results for the decay width are given in the Table 1 for two sets of model parameters. The main 
contribution comes from the box diagram. The contribution from vector mesons has a constructive interference while 
the scalar an contribution has a destructive one. The results are in satisfactory agreement with Crystal Ball data 
0.45 ± 0.12 [Hj and the present value 0.57 ± 0.21 given in PDGpTl. 

It is also very instructive to consider the invariant mass distribution. In Figures [21 and 01 the invariant mass 
distribution of the two-photons is shown for the scalar meson contribution, vector mesons contribution, scalar + 
vector mesons and total. In Figure ^ the results of our calculations of the invariant mass distribution are compared 
with the calculation in the chiral unitary approach [2fl |. 



IV. CONCLUSIONS 



Earlier calculations of the process r\ — > ir°j-f in the NJL model do not include the momentum dependence of quark 
loops and pseudoscalar-axial- vector transitions and are in satisfactory agreement with the GAMS experiment. 
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FIG. 2: Invariant mass distribution of the two-photons of the scalar meson contribution(dots), vector meson contributions(short 
dash), scalar + vector mesons(dash-dot), quark box(long dash) and total(continuous line) for the set I. 




FIG. 3: Invariant mass distribution of the two-photons of the scalar meson contribution(dots), vector meson contributions(short 
dash), scalar + vector mesons(dash-dot) , quark box(long dash) and total(continuous line) for the set II. 



Recently, the new experimental data on this decay have been obtained and the value of the decay width is almost 
two times smaller. A number of theoretical estimates is also obtained, and it seems that the momentum dependence 
of amplitudes is important for a correct description of this process ( "all-order" estimate in ChPT) . 

In the present work, the contributions from quark box, scalar and vector pole diagrams are considered with the full 
momentum dependence. The pseudoscalar-axial-vector transitions are also taken into account. 

The obtained result is consistent with recent experiments, theoretical estimates of ChPT 0, and the chiral 
unitary approach ppj. 

In future, we plan to consider the polarizability of pions and also decays of vector mesons p(oj) — > r](ir)ir~f. 
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FIG. 4: Invariant mass distribution of the two-photons of the total contributions for the set I(dashes), set Il(dots) together 
with the results of the chiral unitary approach |20(| . 



Appendixes 
A. Lagrangian 

Lagrangian Q can be rewritten in the form (see [HI EH) 

1 9 

C = q(id ~ m°)q + - £[gH (q\\q) 2 + G { + ] (q ll5 \\q) 2 

G\ 



where 



-GU (QKq)(q^ s q) + G^ (qi-f 5 \ u q)(qi"/ 5 \ s q) + ^[(^A^) 2 + {qisl ^Xq) 2 ] 



i=0 



X'i = \ (i=l,...,7), A' 8 = A„ = (V2X + Aa)/V3, 

Ag = A s = (-A + V2X 8 )/V3, 

G[ ±] = G { 2 ±] = G { 3 ±] =G± iKm s h{m s ), 

G^ - G^ } = Gf ] = G { ± ] =G± 4Km u h{m u ), 

Gi ±J = G T 4Km s h(m s ), = G, G<±> = ±AV2Km u h(i 



(12) 



(13) 



B. Polarization loops 

Polarization loops in different channels after the PV regularization 

e -izm imj R ..^ = e -izm imj [i _ (i + izA 2 )e' tzA2 



(14) 



take the form (see [22| for the expressions for the polarization loops with equal indices) 

-~ + t^ 2 (! - y 2 ) - 2 [i m i - m if - y( m i - 



n^V) 



n^V) 



n«kV) 
n» 2 ) 



n^( P 2 ), 



„2 



n^ 2 ) = n l #( P 2 )- P 2 y ^(l-y 2 ) y "Raided 



" ^n^( P 2 ) + ^fn^( P 2 ), 



]\T /•! /■<» ,J 7 

n^(p 2 ) + 2m I m J i -| y dy y Hj^.f*)^, 



-R lj (z)e 



A = — (1 - y 2 ) - - [(m t - mj) - y(ra\ - m 2 )] . 



Vertex functions 



The most simple form have the vertex functions for the vector p and the isovector scalar meson do, namely 

V aa = g ao la , V p = gp'i^p^. 
The matrices G and II for dp and p mesons have the form 



G[-\n ao ( P 2 ) = uTs(,P 2 ), 

uu ( r Jl\ 



G p = G 2 ,Tl p (j/)=U^ v (p 2 ) 

For the pion and kaon, additional axial- vector components appear in the vertex function due to pseudoscalar 
vector mixing 



V-k = flW*75(l + AtvP)k, V k = 5^*75(1 + &kp)K 



Here G and II are 



G K = 



l) • n ' (p2) 



n^(p 2 ) n^ u A ( P 2 ) 







j u /n^ P (p 2 ) n^( P 2 ) 



S gJ ,IW) ln^ P (p 2 ) rr^ 2 ), 



5 We suppress flavor indices. 
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Therefore, the vertex function of the 77 meson have four components: strange and non-strange pseudoscalar and 
axial- vector 

Vr, = 5^75(1 + A Vu p)r] u + 5^75(1 + \ s p)Vs = (20) 
= 5^75 (cos Ot,^ - sin6^?7 s + A v p(cosQ v r} u - sin©^)), 

where ©^ and ©^ are the mixing angles for pseudoscalar and axial- vector components. The matrices G and H(p 2 ) 
are four-by-four matrices 

G = (T G 2 ) • G<+> = {oh gF>) ° 2 = di ^(a,&( (21) 



P,A 



D. Amplitudes r\ — ► 77, p — > ??(7r)7 

The amplitude for the two-photon decay width of the pseudoscalar meson has the form 

A(P - 77) = e 2 3 p 77 (M 2 , ql q 2 ) e^ af3 tf<fi , (22) 
where <?i, 52 are the momentum of photons and e 1 , e 2 , are the polarization vectors of the photons, 

g^(M 2 ,q 2 ,q 2 ) = I u (M 2 7 q 2 7 q 2 )g 7:7 (23) 
g vll (M 2 ,q 2 ,q 2 ) = puiM^qlq^g^-^IsiM^qlq 2 ^. 

The loop integrals Ij(Mp) are given by 

1 l-Xl 

7 ^ M -^ 2 ) ^ / _ xi[1 _ xi _ X2)q2 Zlil - Xl - X2 )ql X1X2 M 2 ■ (24) 



The amplitudes for the processes p(ui) — > 77(77)7 have the form 

A(PT/ 7 ) - 5p e 9Pp7 (M|, g 2 , g 2 ) e „„ Q/J g£ , (25) 

here <7i and are the momentum and the polarization vector of p{oj) meson. 

g^ P - ( (M 2 , ql q 2 ) = I U (M 2 , q\, q 2 )g n , g vpl {M 2 ,q 2 ,q 2 ) = 3I u {M 2 ,q 2 ,q 2 ) griu , 

g^ 7 (M 2 , q 2 , q 2 ) = 3I u (M 2 ,q 2 , q 2 )g w , g v ^(M 2 , q\,q 2 ) = I u {M 2 ,q 2 ,q 2 ) gr , u . (26) 
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